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Abstract. We have observed NGC 3603, the most mas- 
sive visible Hff region known in the Galaxy, with 
ANTU(VLTl)/fSAAC in the near-infrared (NIR) J s , H, 
and -Kg-bands. Our observations are the most sensitive 
observations made to date of this dense starburst region, 
allowing us to investigate with unprecedented quality its 
low-mass stellar population. Our mass limit to stars de- 
tected in all three bands is 0.1M Q for a pre-main sequence 
star of age 0.7 Myr. The overall age of the pre-main se- 
quence stars in the core region of NGC 3603 has been 
derived from isochrone fitting in the colour-magnitude di- 
agram, leading to 0.3—1.0 Myr. The NIR luminosity func- 
tions show that the cluster is populated in low-mass stars 
at least down to 0.1M Q . Our observations clearly show 
that sub-solar mass stars do form in massive starbursts. 

Key words: Telescopes: VLT - H-R and C-M diagrams - 
Stars: pre-main sequence - HII regions: NGC3603 - Galax- 
ies: starbursts 



1. Introduction 

NGC 3603 is located in the Carina spiral arm (RA = 11 , 
DEC = -61°) at a distance of - 6 - 7 kpc (De Pree et 
al. 1999, and references therein). It is the only massive, 
Galactic HII region whose ionizing central cluster can 
be studied at optical wavelengths due to only moderate 
(mainly foreground) extinction of Ay ~ 4.5 m (Eisenhauer 
et al. 1998). The OB stars (> 10M Q ) contribute more 
than 2OOOM to the cluster mass. In comparison to the 
Orion Trapezium system, NGC 3603 with its more than 
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50 O and WR stars (Moffat, Drissen & Shara 1994) pro- 
ducing a Lyman continuum flux of 10 51 s _1 (Kennicutt 
1984; Drissen et al. 1995) has about 100 times the ion- 
izing power of the Trapezium cluster. With a bolometric 
luminosity L hol > 1O 7 L , NGC 3603 has about 10% of 
the luminosity of 30 Doradus and looks in many respects 
very similar to its stellar core R136 (Brandl et al. 1996). 
In fact it has been called a Galactic clone of R136 with- 
out the massive surrounding cluster halo (Moffat, Drissen 
k Shara 1994). In many ways NGC 3603 and R136 can 
be regarded as representative building blocks of more dis- 
tant and luminous starburst galaxies (Brandl, Brandner 
& Zinnecker 1999, and references therein). 

Although NGC 3603 has been the target of many re- 
cent ground and space based studies (Eisenhauer et al. 
1998; Brandner et al. 2000) little is known about its con- 
tent of sub-solar mass stars. Our general aim is to investi- 
gate low-mass star formation in the violent environments 
of starburst regions (Brandl, Brandner & Zinnecker 1999). 
Several fundamental questions arise in this context: Does 
the slope of the IMF vary on small scales? Do low-mass 
stars in a starburst event form together with the most 
massive stars or do they form at different times or on 
different timescales? And finally, one might even ask if 
low-mass stars form at all in such environments. 

2. Observations and Data Reduction 

NGC 3603 was observed in the J s = 1.16 — 1.32/im, 
H = 1.50 - 1.80^m, and K s = 2.03 - 2.30^m broad- 
band filters using the NIR camera ISAAC mounted on 
ANTU, ESO's first VLT. The observations were made dur- 
ing the 4 nights of April 4-6 and 9, 1999, in service mode 
when the optical seeing was equal or better than 0."4 on 
Paranal. Such seeing was essential for accurate photom- 
etry in the crowded cluster and increased our sensitivity 
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Fig. 1. Three colour image of NGC 3603 composed from J s (blue), H (green), and if s -band (red) images. Intensities 
are scaled in logarithmic units; FOV is 3.'4 x 3.'4 (6.2 x 6.2 parsec 2 ). North is up, East to the left. The insert to the 
lower right is a blow up of the central parsec 2 . 



to the faintest stars. The majority of our data were taken 
under photometric conditions. 

Our observing strategy was to use the shortest possible 
frame times (1.77s) to keep the number of saturated stars 
to a minimum. However, due to the system's excellent 
sensitivity, about two dozen of the brightest stars ended 
up being saturated. Nevertheless, this does not impose 
a problem to our study of the low mass stars. Thirty- 



four short exposures were co-added to an effective one- 
minute exposure, the minimum time per pointing required 
to stabilize the telescope's active optics control system. 
Between the 1 minute pointings we moved the telescope 
by up to 20" offsets in a random pattern. This approach 
has several advantages: 

— Enlargement of the observed field of view (FOV) with 
maximum signal-to-noise (S/N) in the cluster center. 
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Fig. 2. a-c. J s versus J s — K s colour-magnitude diagrams of NGC 3603. a contains all stars detected in all three 
wavebands within the central r < 33" (lpc) [1 /D ], b shows the field stars at r > 75" (2.25 pc) [6.3' n ] around the 
cluster, and c shows the cluster population within r < 33" with the field stars statistically subtracted. The dashed 
horizontal line indicates the detection limit of the previous most sensitive NIR study by Eisenhauer et al. (1998). c 
also shows the theoretical isochrones of pre-main sequence stars of different ages from Palla & Stahler (1999) and the 
main sequence for dwarfs. For comparison we've plotted some corresponding stellar masses next to the isochrones. 



— Reduction of residual images and other array artefacts, 
using the median filtering technique. 

— Derivation of the "sky" from the target exposures us- 
ing the median filtering technique. No additional time 
for "blank" sky frames outside the cluster was re- 
quired. 

The sky frames have been computed using between 15 
and 37 subsequent exposures per waveband and night, 
and careful eye-inspection showed that all sources have 
been efficiently removed using our modified median filter- 
ing technique which returns the lower 1/3 instead of the 
mean (1/2) value. We subtracted the sky-background and 
flat-fielded each exposure using the twilight flat-fields pro- 
vided by ESO. The relative position offsets were derived 
from cross-correlating the images; the exposures were co- 
aligned on a 0.5 x 0.5-pixel subgrid for better spatial reso- 
lution, and then added together using the median filtering 
technique. The resulting images are 3.'4 x 3.'4 in size with 
pixels of 0."074. The effective exposure times of the final 
broadband images in the central 2.'5 x 2.'5 is 37, 45, and 48 
minutes in J s , H, and K s , respectively. Fig. |l] shows the 
impressive 3-colour composite image. The brightest star 
in the FOV (80" northeast of the core) is the red super- 
giant IRS 4 (Frogel, Persson, & Aaronson 1977). The ring 
nebula and the bipolar outflows around the blue super- 
giant Sher 25 (Brandner et al. 1997) about 18" north of 



the core are clearly detected in the K band. The image 
also shows the three proplyd-like objects that have been 
recently discovered by Brandner et al. (2000); these pro- 
plyds are similar to those seen in Orion but about 20-30 
times more extended. About 1' south of the central cluster, 
we detect the brightest members of the deeply embedded 
proto cluster IRS 9. 

In order to derive the photometric fluxes of the stars 
we used the IRAF implementation of DAOPHOT (Stet- 
son 1987). We first ran DAOFIND to detect the individual 
sources, leading to ~20,000 peaks in each waveband. Many 
of these may be noise or peaks in the nebular background 
and appear only in one waveband. In order to reject spu- 
rious sources, we required that sources be detected inde- 
pendently in all three wavebands, and that the maximal 
deviation of the source position centroid between different 
wavebands be less than 0.075". The resulting source list 
contains 6967 objects in the entire FOV. 

We then flux-calibrated the images using the faint NIR 
standard stars from the lists by Hunt et al. (1998) and 
Persson et al. (1998). Because of the stringent require- 
ments on the seeing the PSF did not noticeably change 
during our observations and the systematic photometric 
errors are dominated by uncertainties in the aperture off- 
sets. (A detailed error analysis will be part of a subsequent 
paper). Comparing our photometric fluxes of numerous 
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Fig. 3. The observed luminosity functions of stars de- 
tected at all 3 wavebands in the central stellar cluster with 
r < 33" (lpc). The magnitudes are apparent J s (solid), 
H (dotted), and K s (dashed) magnitudes. The number 
counts have not been corrected for incompleteness. 



sources with the fluxes derived by Eisenhauer et al. (1998) 
yields a systematic offset of 0.1 m in J s and 0.05 m in K s . 



3. Results and Interpretation 

Fig. 0a-c shows the resulting colour-magnitude diagrams 
(CMD). The left plot contains all stars detected in all 3 
wavebands within r < 33" (1 pc). Since NGC 3603 is 
located in the Galactic Plane we expect a significant con- 
tamination from field stars. To reduce this contribution 
we followed a statistical approach by subtracting the av- 
erage number of field stars found in the regions around 
the cluster at r > 75" (central plot) per magnitude and 
per colour bin (0.5 mag each). The accuracy of our statis- 
tical subtraction is mainly limited by three factors: first, 
we cannot rule out that low-mass pre-main sequence stars 
are also present in the outskirts of the cluster. Second, be- 
cause of crowding on one hand and dithering which leads 
to shorter effective integration times outside the central 
2.'5 x 2.'5 on the other hand, the completeness limit varies 
across the FOV. Third, local nebulosities may hide back- 
ground field stars. However, none of these potential errors 
affects our conclusions drawn from the CMD. 

The resulting net CMD for cluster stars within r < 33" 
of NGC 3603 is shown at the right in Fig. g. We overlayed 
the theoretical isochrones of pre-main sequence stars from 
Palla & Stahler (1999) down to 0.1Af Q . We assumed a dis- 
tance modulus of (m — M) = 13.9 based on the distance 
of 6 kpc (De Pree et al. 1999) and an average foreground 
extinction of Ay — 4.5 m following the reddening law by 
Rieke & Lebofski (1985). Applying the isochrones to mea- 
sured magnitudes could be misleading since the theoret- 
ical calculations include only stellar photospheres while 
the stars still may be surrounded by dust envelopes and 
accretion disks. This would lead to excess emission in the 



NIR and make the stars appear younger than they actu- 
ally are. Typical excess emission of classical T Tauri stars 
in the Taurus-Auriga complex have been determined as 
AH = 0.2 m , and AK = 0.5 m (Meyer, Calvet & Hillen- 
brand 1997). The upper part of the cluster- minus-field 
CMD clearly shows a main sequence with a marked knee 
indicating the transition to pre-main sequence stars. The 
turn-on occurs at J s « 15. 5 m (to « 2.9Af Q ). Below the 
turn-on the main-sequence basically disappears. We note 
that the width of the pre-main sequence in the right part of 
Fig. |2j does not significantly broaden toward fainter mag- 
nitudes, indicating that our photometry is not limited by 
photometric errors. The scatter may in fact be real and 
due to varying foreground extinction, infrared excess and 
evolutionary stage. In that case the left rim of the dis- 
tribution would be representative of the "true" colour of 
the most evolved stars while the horizontal scatter would 
be primarily caused by accretion disks of different incli- 
nations and ages. Fitting isochrones to the left rim in the 
CMD yields an age of only 0.3 — 1.0 Myr. Our result is 
in good agreement with the study by Eisenhauer et al. 
(1998) but extends the investigated mass range by about 
one order of magnitude toward smaller masses. Because 
of ~ 10 magnitudes range in luminosities in the crowded 
core region our sensitivity limits of J s w 21 m , H ss 20 m 
and K s « 19 m don't appear to be exceedingly faint (and 
are about 3 magnitudes above ISAAC'S detection limit 
for isolated sources). However, only VLT/ISAAC's high 
angular resolution, PSF stability, and overall sensitivity 
enabled us to study the sub-solar stellar population in a 
star burst region on a star by star basis. 

Fig. shows the luminosity functions for stars de- 
tected in all 3 wavebands. For the purpose of this letter 
we have not attempted to correct our number counts for 
incompleteness, i.e., an increasingly significant percent- 
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age of stars will be undetected toward fainter magnitudes. 
Thus we have no means to determine whether the appar- 
ent turnover at K a f=a 16. 5 m is real or an observational 
artifact, but we do note that the mass spectrum is well 
populated down to K s ~ 19 m , corresponding to a 0.1M Q 
pre-main sequence star of 0.7 Myr (Zinnecker, McCaugh- 
rean & Wilking 1993). A detailed analysis of the low-mass 
IMF will be the subject of a subsequent paper. 
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